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INTRODUCTION 

The  objective  of  chis  study  was  to  evaluace  the  ability  of  synthetic 
aperture  radar  (SAR)  data  collected  by  the  SEASAT  satellite  to  provide 
information  on  bottom  topography  in  oceanic  areas.  Investigations  were 
undertaken  to  evaluate  SEASAT 's  ability  to  detect  uncharted  or  mispositioned 
submerged  features  hazardous  to  marine  navigation,  to  position  these 
features  accurately,  and  to  investigate  whether  or  not  depth  Information 
can  be  extracted  from  SAR  data.  Specifically,  three  analyses  were  per¬ 
formed.  First,  correlations  were  made  between  surface  patterns  present  on 
SAR  imagery  in  oceanic  areas  and  hydrographic  features.  Secondly,  tech¬ 
niques  to  enhance  subtle  SAR  signatures  Indicative  of  bottom  features  were 
evaluated.  And  finally,  several  techniques  to  measure  relative  radar  back- 
scatter  return  were  evaluated  in  an  attempt  to  better  quantify  the  relation¬ 
ship  between  the  surface  features  and  the  related  bottom  topography. 

The  data  studied  in  this  report  was  collected  by  the  SEASAT 
satellite.  Among  the  instrumentation  carried  by  SEASAT,  which  was 
launched  during  June  of  1978,  was  a  synthetic  aperture  radar  (SAR). 

This  satellite  collected  over  500  passes  of  SAR  data  before  suffering 
a  catastrophic  power  loss  in  October  of  1978.  The  SAR  on  board  SEASAT 
was  an  L-band  (23.5  cm  wavelength)  radar.  It  collected  25  x  25  meter 
resolution  imagery  with  a  ground  swath-width  of  100  km  and  a  length 
of  up  to  4000  km. 

The  wavelength  at  which  the  SEASAT  SAR  system  operates  (23.5  cm) 
does  not  penetrate  more  than  a  few  centimeters  into  the  water  surface. 
Therefore,  the  recorded  return  signal  is  primarily  reflected  from  the 
water's  surface  and  is  a  function  of  the  surface  roughness  and  the 
surface  slope.  The  surface  features  reported  in  this  study  are  believed 
to  be  the  result  of  a  hydrodynamic  interaction  between  a  physical 
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oceanic  process  and  a  bottom  topographic  feature.  This  interaction 
results  in  an  alteration  of  either  the  surface  roughness  or  slope; 
hence  a  surface  anomaly  is  observed  on  the  SAR  imagery. 

Five  different  areas  from  four  SEASAT  orbits  were  studied. 

Initial  comparisons  between  the  surface  features  present  in  the 
coastal  regions  of  these  passes  and  bottom  topography  information  were 
made  to  see  if  any  correlation  existed.  In  a  very  high  percentage  of 
the  cases,  a  correlation  did  exist  between  a  surface  signature  on  the 
SAR  imagery  and  a  bottom  topographic  feature.  During  this  analysis  it 
was  observed  that  in  some  areas  the  SAR  surface  signatures  were 
distinct,  while  in  other  areas  the  surface  features  were  quite  subtle. 

In  those  areas  of  subtle  surface  features,  a  second  analysis  was 
conducted  to  see  if  these  subtle  features  could  be  made  more  visible 
through  photographic  enhancement  and  SAR-processing  techniques.  It  was 
thought  that  these  special  processing  techniques  might  also  bring 
out  features  which  were  not  previously  visible  on  the  survey  processed 
(quick-look)  imagery.  The  photographic  techniques  included  photographic 
enlargements  and  use  of  different  levels  of  contrast  on  the  photographic 
paper.  The  SAR  processing  techniques  involved  special  optical  processing 
of  the  SAR  signal  film  including  motion  corrections,  looking  at  the 
proper  dynamic  range  to  display  oceanographic  applications,  and 
variation  of  recording  film  types  to  test  several  film  contrasts. 

Relative  radar  backscatter  measurements  were  made  in  the  third 
analysis  in  an  attempt  to  see  if  a  quantitative  relationship  existed 
between  the  observed  SAR  radar  returns  and  water  depth.  These  rela¬ 
tive  backscatter  measurements  were  obtained  from  the  radar  signal 
film  using  the  ERIM  hybrid  digitizer  and  from  JPL  digitally  processed 
data. 
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SUMMARY 


This  study  evaluated  the  use  of  synthetic  aperture  radar  (SAR) 
imagery  collected  by  the  SEASAT  satellite  to  provide  bathymetric  in¬ 
formation.  The  SAR  system  can  provide  bathymetric  information  by 
imaging  surface  perturbations  which  are  a  result  from  a  hydro- 
dynamic  interaction  between  a  physical  oceanic  process  (such  as 
a  current  or  wave)  and  a  distinct  hydrographic  feature. 

Imagery  from  four  SEASAT  passes  was  used  to  study  five  study 
sites.  These  sites  included:  the  DMA-Bahama  Photobathymetric 
Calibration  Area  in  the  Northwestern  Bahamas;  a  portion  of  the  South¬ 
western  Great  Bahama  Banks,  just  north  of  Cuba;  the  Nantucket  Shoals 
region  of  Cape  Code;  portions  of  the  English  Channel;  and  the  area 
surrounding  North  Rona  Rock,  off  the  coast  of  Scotland. 

The  surface  features  present  in  the  two  test  sites  in  the  Bahamas 
showed  a  high  degree  of  correlation  to  bottom  topographic  features 
(60  to  100%),  but  were  quite  subtle,  and  were  located  only  at  the  edges 
of  the  Great  and  Little  Bahama  Banks.  Special  photographic  and  optical 
processing  techniques  were  used  on  these  data  to  find  other  depth- 
related  surface  patterns  but  did  not  bring  forth  any  additional  infor¬ 
mation. 

The  surface  patterns  present  in  the  other  study  sites  were  quite 
distinct  and  an  examination  of  the  position  of  these  patterns 
revealed  a  strong  correlation  (77%  or  higher)  relative  to  the 
position  of  distinct  bottom  topographic  features,  as  obtained  from 
bathymetric  charts. 

Relative  radar  backscatter  measurements  were  obtained  at  four 
of  the  five  test  sites  (the  Southern  3ahamas  Test  Site  was  excluded). 
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These  backscatter  measurements  were  produced  in  transect  form  and 
were  compared  to  bathymetric  cross-section  transects  of  Che  same 
areas.  A  high  percentage  of  the  major  bottom  topographic  features 
on  the  bathymetric  cross-sections  were  positioned  directly  below  a 
major  surface  feature  on  the  SAR  data  (which  is  indicated  by  a 
significant  (3-6  dB]  drop  or  rise  on  the  radar  backscatter  transect). 

Also  presented  is  ancillary  information  (such  as  tidal  and 
meteorological  data)  and  Landsat  imagery  which  support  the  observa¬ 
tion  chat  the  surface  patterns  present  on  the  SAR  imagery  are  an 
indirect  result  of  a  hydrodynamic  interaction  between  a  physical 
oceanic  process  (i.e.,  a  current)  and  a  distinct  bottom  topographic 
feature. 

This  study  clearly  demonstrates  the  feasibility  of  utilizing 
synthetic  aperture  radar  imagery  to  provide  bathymetric  information. 

SAR  can  be  used  to  detect  the  presence  of  certain  types  of  bathymetric 
features  and  can  provide  a  rough  indication  of  their  location.  However, 
it  should  be  noted,  the  limitations  of  SAR  to  obtain  this  information 
(i.e.,  under  what  oceanic  and  meteorological  conditions)  is  not  totally 
understood  at  present. 
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BACKGROUND 

During  examination  of  SEASAT  Synthetic  Aperture  Radar  (SAR) 
imagery,  it  was  noted  that  a  variety  of  surface  anomalies  existed, 
both  in  deep  and  shallow  water  areas.  A  theory  has  been  postulated 
that  many  of  these  anomalies  are  a  result  of  an  interaction  between  a 
physical  oceanic  process  and  a  distinct  bottom  feature,  and  as  such, 
represent  a  way  of  using  SAR  imagery  to  gain  bathymetric  information. 

Five  different  physical  oceanographic  phenomena  have  been 
observed  on  SAR  imagery  that  permit  detection  of  underwater  hazards. 

These  include: 

1.  The  change  in  direction  and  wavelength  of  ocean  swell  as  they 
enter  coastal  regions  which  can  be  observed  on  SAR  imagery; 

2.  Observation  of  classic  internal  wave  patterns  over  continen¬ 
tal  shelf  regions; 

3.  Anomalous  SAR  backscatter  signatures  found  in  deep  ocean 
regions  over  topographic  "bumps"  on  the  ocean  bottom, 

4.  A  distinct  change  in  radar  backscatter  which  is  associated 
with  an  ocean  swell  propagating  over  a  distinct  depth  dis¬ 
continuity;  and 

5.  A  distinct  change  in  radar  backscatter  which  is  associated 
with  a  current  flowing  over  a  bottom  feature. 

Models  which  incorporate  the  way  an  ocean  waves  changes  as  it 
enters  shallow  waters  can  be  used  to  gain  an  estimate  of  water  depths 
in  coastal  regions.  When  a  gravity  wave  enters  a  shallow-water  region 
from  deep  water,  its  wavelength  and  direction  change  in  direct  proportion 
to  the  water  depth.  By  measuring  the  change  in  wavelength  or  direction, 
an  estimate  of  the  water  depth  can  be  obtained.  Since  a  SAR  can  measure 
both  ocean  wavelength  and  direction,  this  information  can  be  used  in  a  wave 
refraction  model.  This  was  done  in  studies  by  Shuehman,  et  al.  (1979). 
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Basic  SAR/Oceanographic  theory  holds  that  the  microwave  energy 
from  the  radar  is  being  reflected  from  the  ocean's  surface  by  the  small 
capillary  and  ultra-gravity  waves  present  at  the  ocean's  surface. 

Capillary  and  ultra-gravity  waves  have  wavelengths  between  one  and  fifty 
centimeters  and  are  generated  by  winds.  Any  change  in  this  small  wave 
structure  will  alter  the  radar  backscatter  being  received  by  the  SAR. 

Internal  wave  packets  propagating  shorewards  over  a  continental 
shelf  region  act  as  a  mechanism  permitting  SAR  detection  of  changes  in 
bottom  form.  The  internal  wave  energy  field  alters  the  small  capillary 
and  ultra-gravity  waves  sufficiently  so  that  the  internal  waves  are 
visible  on  SAR  imagery  (Shuchman  and  Kasischke,  1979).  An  example  of 
this  phenomena  is  presented  in  Figures  1  and  2.  Surface  anomalies,  similar 
to  internal  waves,  but  not  as  well  structured  (into  wave  packets),  have 
been  noticed  on  SAR  images  in  deep-water  areas  between  major  continental 
shelf  boundaries  (Shuchman  and  Kasischke,  1979).  The  interesting  fact 
about  these  surface  anomalies  is  that  in  large  part,  they  occurred  over  a 
major  underwater  ridge  or  seamount  (see  Figures  3  and  4). 

DeLoor  and  VanHulten  (1978)  observed  bottom-related  surface  features 
on  radar  imagery  collected  over  shallow  water  regions  of  the  North  Sea. 

3.1  SAR  CONCEPTS 

A  SAR  system  offers  several  advantages  over  passive  remote 
sensing  techniques.  Synthetic  Aperture  Radar  is  an  active  imaging 
device  that  senses  the  environment  with  short  electromagnetic  waves. 

As  active  sensors,  radars  provide  their  own  illumination  in  the  micro¬ 
wave  region  of  the  electromagnetic  spectrum  and  thus  are  not  affected 
by  diurnal  changes  in  emitted  or  reflected  radiation  from  the  earth's 
surface.  Additionally,  SAR's  have  the  recognized  advantage  of  being 
able  to  image  the  earth's  surface  independent  of  weather  conditions  or 
solar  illumination  and  to  provide  synoptic  views  of  the  ocean  at  high 
resolution.  Most  radars  operate  in  the  frequency  region  of  300  MHz 
(1  m)  to  30  GHz  (1  cm),  and  bandwidths  within  this  region  are  commonly 
designated  by  letters.  The  SEASAT  SAR  discussed  in  this  paper  is 
designated  as  an  L-Band  (23.5  cm)  system. 
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1  Lit  ion  1232,  21  September  1978. 


Figure  3.  SAR  Imagery  of  Surface  Anomolies  Present  During 
SEASAT  Revolution  762,  19  August  1978,  Between 
Scotland  and  Iceland. 
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Synthetic  Aperture  Radar  is  a  coherent  airborne  or  spaceborne 
radar  that  uses  the  motion  of  a  moderately  broad  physical  antenna 
beam  to  synthesize  a  very  narrow  beam  thus  providing  fine  azimuthal 
(along-track)  resolution  (Brown  and  Porcello,  1969;  Harger,  1970). 

Fine  range  (cross-track)  resolution  is  achieved  by  transmitting  either 
very  short  pulses  or  longer  coded  pulses  which  are  compressed  by 
matched-filtering  techniques  into  equivalent  short  pulses. 

Usually,  the  coded  pulse  is  a  waveform  linearly  modulated  in  fre¬ 
quency  . 

Typically,  the  phase  history  of  a  scattering  point  in  the  scene 
is  recorded  on  photographic  film  as  an  anamorphic  (astigmatic)  Fresnel 
zone  plate.  The  parameters  of  the  zone  plate  are  set  in  the  azimuth 
direction  by  the  Doppler  frequencies  produced  by  the  relative  motion 
between  the  sensor  and  the  point  scatterer,  and  in  the  range  direction 
by  the  structure  of  the  transmitted  pulses.  The  film  image  is  a 
collection  of  superimposed  zone  plates  representing  the  collection  of 
point  scatterers  in  the  scene.  This  film  is  used  by  a  coherent  optical 
processor  which  focuses  the  anamorphic  zone  plates  into  the  points 
which  produced  the  microwave  scatter  of  the  scene  (Kozma,  et  al.,  1972). 
Recently,  SAR  processing  has  employed  digital  techniques. 

However,  when  using  a  SAR  system  to  image  moving  targets  such  as 
the  ocean  surface,  unique  problems  occur  when  correlating  the  signal 
data.  Because  moving  targets  (such  as  waves  or  currents)  perturb  the 
Doppler  frequencies,  and  hence  the  phase  histories  recorded  by  the 
signal  receiver,  conventional  processing  of  these  signal  histories 
produces  images  of  the  ocean  that  are  defocused  relative  to  a  stationary 
target.  Defocusing  by  the  along-track  (azimuth)  velocity  of  moving 
ocean  surfaces  can  be  refocused  by  re-adjusting  the  azimuth  focus 
(cylindrical  focus)  an  amount  proportional  to  the  relative  velocity 
of  the  ocean  surface. 
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Similarly,  the  radial  morion  of  a  moving  ocean  surface  imaged  by 
a  SAR  will  also  perturb  the  signal  history  of  a  scatterer.  Radial 
velocity  (motion  towards  the  radar  look,  direction)  produces  an 
apparent  tilt  to  the  phase  history  as  recorded  by  the  signal  receiver. 
Additionally,  the  scatterer  history  will  also  shift  across  the  signal 
record.  This  is  referred  to  as  "range  walk"  and  can  be  compensated 
for  by  a  rotation  of  the  cylindrical  optics  in  the  processing  of  the 
signal  histories. 

The  signal  histories  recorded  by  a  SAR  frequently  have  a  range  of 
return  brightness  spanning  60  decibels  (dB).  This  range  of  signal  returns 
exists  whether  the  histories  are  recorded  in  a  digital  or  optical  format 
(on  SAR  signal  film) .  Problems  arise  when  converting  the  signal  histories 
to  ouptut  images  in  either  an  optical  or  digital  format. 

If  the  output  product  desired  is  an  image  film,  it  must  be  realized 
that  most  photographic  emulsions  can  only  record  light  intensities  with  a 
dynamic  range  of  17-20  dB  with  a  linear  response.  Further,  because  of 
brightness  compressions  associated  with  this  non-linear  response  pattern,  many 
small  reflectivity  variations  in  an  image  can  be  lost  due  to  "lumping"  of 
brightness  levels  at  the  extreme  ends  of  the  intensity  response  window. 

Thus  a  film  representation  of  SAR  output  imagery  can  "mash"  features  detected 
by  the  SAR. 

The  large  dynamic  range  scenes  imaged  by  a  SAR  can  be  displayed  on 
film  bv  creating  a  number  of  output  films  encompassing  the  entire  dynamic 
range  of  a  scene  imaged  by  SAR  (Larson,  et  al.,  1975).  If  computer 
compatible  tapes  (CCT's)  are  used  to  store  SAR  output  then  approximately 
40  dB  worth  of  information  can  be  retained. 

3.2  SEASAT  SAR  DATA 

3efore  entering  into  a  discussion  of  the  various  measurement  and 
analysis  techniques  used  during  this  project,  we  should  first  mention 
the  various  forms  the  SEASAT  SAR  data  are  available  to  the  user. 

SEASAT  SAR  data  can  be  obtained  in  three  ways:  1)  JPL  survey  pro¬ 
cessed  imagery;  2)  SAR  Doppler  signal  histories:  and,  3)  digital 
imagery  and  capes. 
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JPL  survey  processed  data  is  available  Co  the  user  in  either 
a  master  positive  or  a  master  negative  fora.  Each  pass  of  SEASAT 
data  covers  four  sections  of  70  mm  wide  film,  with  each  section  con¬ 
taining  a  25  kilometer  swath  of  data.  The  scale  factor  of  this  imagery 
is  approximately  1:500,000  with  a  resolution  of  30-40  meters.  This 
data  is  relatively  inexpensive  and  can  usually  be  obtained  in  a  few 
weeks.  This  data  has  been  nominally  geometrically  corrected  and 
usually  is  adequate  for  a  quick  first-look.  The  disadvantages 
associated  with  this  data  type  are:  1)  the  imagery  is  a  third 
generation  film  product;  2)  the  data  is  not  always  optimally  pro¬ 
cessed  with  respect  to  the  proper  dynamic  range  for  oceanographic 
observations  and  the  full  dynamic  range  of  the  data  is  unexploited; 

3)  a  limited  number  of  the  total  SEASAT  passes  are  available  in  this 
form;  4)  the  imagery  has  not  been  processed  to  account  for  ocean 
surface  motion;  and,  5)  the  full  potential  resolution  of  the  SAR  system 
is  not  realized. 

The  second  type  of  SEASAT  SAR  data  utilized  in  this  study  was  the 
SAR  Doppler  signal  histories  recorded  on  film.  The  difference  between 
signal  film  and  image  film  should  be  noted  at  this  time  to  avoid  con¬ 
fusion.  Signal  film  is  used  to  record  the  phase  history  of  the  various 
scattering  points  on  the  ground  being  illuminated  by  the  SAR.  By 
passing  a  collimated  beam  of  monochromatic  light  through  this  signal  film 
within  an  optical  processing  system,  the  ground  image  is  reconstructed. 

By  placing  a  recording  film  in  the  output  plane,  an  image  film  is 
created . 

The  advantages  of  using  SAR  signal  film  lies  in  the  fact  that 
one  can  custom  process  the  data  according  to  the  needs  of  the  analyst. 

The  data  can  be  geometrically  rectified,  the  proper  dynamic  range  of 
the  data  for  oceanographic  applications  can  be  used,  maximum  resolution 
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can  be  obtained  (10  meters) ,  and  any  image  films  are  first  genera¬ 
tion  film  products.  In  addition,  other  output  products  are  avail¬ 
able.  Digitized  data  (on  CCT's)  can  be  obtained  and  relative  radar  back- 
scatter  measurements  can  be  made  by  placing  a  recording  densitometer 
in  the  output  plane  of  the  optical  processor.  The  main  disadvantages 
of  this  data  source  are  that  special  sophisticated  processing  facilities 
are  needed,  and  associated  with  these  facilities  are  considerable  man¬ 
power  expenses . 

The  last  source  of  SEASAT  SAR  data  available  are  digitally  pro¬ 
cessed  computer  compatible  capes.  These  tapes  are  typically  available 
from  the  Jet  Propulsion  Laboratory  ( JPL) .  Each  JPL  CCT  contains  a  single 
100  x  100  kilometer  scene  and  the  25  meter  resolution  data  has  not  been 
geometrically  rectified.  The  advantages  of  this  data  source  is  that 
close  to  the  full  (up  to  37  dB)  dynamic  range  of  the  data  can  be  ex¬ 
ploited,  very  high  quality  photo  imagery  can  be  produced,  and  a  variety 
of  computer  based  analyses  can  be  performed.  The  JPL  CCT's  dis¬ 
advantages  are:  1)  the  data  has  not  been  geometrically  rectified; 

2)  the  scene  sizes  are  only  100  x  100  kilometers,  with  a  limited 
number  of  scenes  available;  3)  Che  amount  of  information  (36  million 
bytes)  in  each  scene  requires  a  large  computer,  and  because  of  this 
large  size,  processing  is  quite  expensive;  and,  4)  processing  of  this 
data  requires  a  set  of  sophisticated  software  programs. 
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DATA  SETS 

The  choice  of  test  sites  was  a  compromise  between  1)  selecting 
enough  data  sets  to  adequately  sample  the  variety  of  conditions  which 
produced  (apparent)  bottom  related  surface  signatures;  2)  the  avail¬ 
ability  of  SEASAT  SAR  data  in  coastal  regions;  3)  the  different  tech¬ 
niques  available  to  analyze  the  data,  and  the  different  data  products 
available  of  the  test  sites;  and,  4)  the  availability  of  resources  to 
adequately  analyze  the  data  as  outlined  in  the  statement  of  work. 

Accordingly,  five  test  sites  were  chosen  (see  Figure  5): 

Site  A:  DMA-Bahamas  Photobathymetric  Calibration  Area.  The  area  in 
and  around  Bimini  Island  (the  northwestern  corner  of  the  Great  Bahama  Bank) 
and  the  area  in  and  around  the  western  tip  of  Grand  Bahama  Island  (the 
southwestern  comer  of  the  Little  Bahama  Bank) . 

Site  B:  Southern  Bahamas.  The  southwestern  corner  of  the  Great 
Bahama  Banks . 

Site  C:  Nantucket  Island.  The  shoal  area  immediately  east  of 
Nantucket  Island,  Massachusetts,  including  the  Nantucket  Shoals  and 
the  area  just  southeast  of  Cape  Cod. 

Site  D:  English  Channel.  The  northern  area  of  the  Straits  of 
Dover  and  the  Thames  Estuary  region. 

Site  E:  North  Rona  Rock.  An  area  northwest  of  Cape  Wrath, 

Scotland,  surrounding  the  island  referred  to  as  North  Rona  Rock. 

Data  from  four  different  SEASAT  orbits  were  used  for  the  various 
analyses.  Study  Sites  A  and  B  were  covered  by  two  different  SEASAT 
orbits:  Revolution  407,  an  ascending  orbit,  passed  over  the  two  test 
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sices  at  approximately  12:46  (GMT)  on  25  July  1978;  and  Revolution  651 
was  collected  on  11  August  1978  at  approximately  12:26  (GMT).  Study 
Site  C  was  covered  during  SEASAT  Revolution  880  (an  ascending  pass)  on 
27  August  1978  at  12:34  (GMT).  Finally,  data  for  Study  Sites  D  and  E 
were  collected  during  SEASAT  Revolution  762  (an  ascending  pass)  on 
19  August  1978  at  6:47  (GMT). 

Table  1  gives  a  summary  of  the  Study  Sites,  the  SEASAT  Revolution(s) 
used  for  each  Site,  the  data  types  available  (for  this  study)  for  each 
SEASAT  Revolution,  and  a  summary  of  the  analysis  techniques  used. 
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TABLE  1.  SUMMARY  OF  DATA  TYPES  AND  DATA  ANALYSIS  TECHNIQUES 
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te  E  -  Nor Lli  Kona  762  Signal  Film  1.  Bathymetric  Correlation  (fror.  image  film) 

goc|t  2.  Hybrid  Digitizer  Backscatter  Measurements 
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BATHYMETRIC  CORRELATIONS 

The  first  analysis  performed  on  each  SAR  data  set  was  to  overlay 
a  clear  transparency  reproduction  of  a  hydrographic  chart  onto  the  SAR 
imagery  to  see  how  closely  the  surface  signatures  correlated  to 
bathymetric  features  (this  process  will  be  referred  to  as  bathymetric 
correlation) .  The  purpose  of  these  bathymetric  correlations  was  to  see 
if  the  observed  features  occurred  near  a  distinct  bottom  feature.  What 
was  looked  for  in  this  analysis  were  general  correlations.  If  a  sur¬ 
face  feature  occurred  near  (i.e.,  within  a  few  hundred  meters),  but  not 
exactly  above  a  bottom  topographic  feature,  a  correlation  was  said  to 
exis  t . 

Figure  6  shows  a  photo  mosaic  made  from  the  four  sub-swaths  of  the 
JPL  produced  imagery  of  SEASAT  Revolution  880,  over  Study  Site  C.  Figure 
7  is  a  blow-up  of  the  area  around  Nantucket  Island  with  a  hydrographic 
chart  overlay.  These  figures  serve  to  illustrate  both  the  type  of  SAR 
imaged  oceanic  surface  features  being  studied,  and  the  degree  of  correla¬ 
tion  which  existed  between  these  features  and  bottom  topographic  infor¬ 
mation. 

Table  2  summarizes  the  bathymetric  correlations.  Given  in  Table 
2  are  the  Study  Site,  SEASAT  Revolution  used  and  percent  correlation. 
Percent  correlation  was  obtained  for  a  given  Study  Site  and  Revolution 
by  comparing  the  positions  of  the  surface  features  on  the  SAR  imagery 
to  the  position  of  distinct  bottom  topographic  features  from  bachymetric 
charts  and  determining  the  ratio  of  the  SAR  features  which  occurred 
over  a  bathymetric  feature.  For  example,  at  Study  Site  C  (Nantucket 
Shoals),  the  77  percent  correlation  figure  was  obtained  by  observing 
that  36  out  of  the  47  distinct  surface  patterns  noted  on  the  SAR 
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Figure  6.  SAR  Imaged  Surface  Features  Near  Nantucket  Island 

(Study  Site  C)  From  SEASAT  Revolution  880,  27  August 
1978.' 
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TABLE  2.  SUMMARY  OF  BATHYMETRIC  CORRELATIONS 


STUDY  SITE 

SEASAT 

REVOLUTION 

PERCENT 

CORRELATION 

A  -  DMA- Bahamas  Photo- 

bacayuietric  Calibration 

Area 

407 

60 

651 

80 

B  -  Southern  Bahamas 

407 

100 

651 

100 

C  -  Nantucket  Shoals 

880 

77 

D  -  English  Channel 

762 

Sub-Site  1:  Thames  River  Estuary 

93 

Sub-Site  2:  English  Channel 

100 

E  -  North  Roaa  Rock 


762 


100 
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imagery  (see  Appendix  A)  occurred  over  a  distinct  bottom  topographic 
feature.  A  complete  discussion  of  this  procedure  and  a  summary  of 
the  individual  data  sets  are  presented  in  Appendix  A. 

The  next  portion  of  this  section  will  present  a  brief  written 
summary  of  bathymetric  correlations  along  with  some  examples  of  the 
SAR  image  features. 

The  bottom-related  surface  features  at  Study  Sice  A  were  along 
the  edges  of  the  Great  Bahama  Bank  (near  Bimini  Island)  and  the  Little 
Bahama  Bank  (near  Grand  Bahama  Island) .  Figure  8  shows  an  example  of 
a  surface  feature  (near  Bimini  Island)which  were  typical  of  this 
area.  No  bottom-related  surface  features  were  noted  on  either 
Bahama  Banks  excepc  near  the  Isaac  Island  group  on  Revolution  651 
imagery,  and  near  Bimini  Island.  The  shallow,  tidal  flat  area  on  the 
southeast  side  of  Bimini  was  evident  on  the  SAR  imagery,  as  were  the 
line  of  scattered  rocks  and  Islets  just  south  of  Bimini. 

The  surface  features  noted  on  SAR  imagery  from  Site  B  (Southern 
Bahama,  see  Figure  91  were  quite  similar  to  those  from  Site  A.  Again, 
surface  features  were  present  near  the  edge  of  the  Great  Bahama  Bank, 
just  north  of  Cuba,  on  both  passes  of  SEASAT  imagery  over  this  area. 

It  should  be  noted  that  the  bottom-related  surface  features  in  Study 
Sices  A  and  B  were  quite  subtle,  which  led  to  the  second  analysis, 
an  attempt  to  enhance  these  subtle  features  through  photographic  and 
optical  processing  techniques.  This  analysis  will  be  presented  in 
the  next  section. 

Almost  8G/o  of  the  suriace  reatures  presenc  on  the  imagery  from 
Sice  C  (Nantucket  Shoals)  were  over  or  near  a  distinct  bottom  feature 
(see  Figures  6  and  7).  This  correlation  improves  to  94%  when  only 
those  features  in  water  less  chan  20  meters  in  depth  were  considered. 
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Bahama  Bank  From  SEASAT  Kevolution  407,  25  July  1978. 


Study  Site  D  was  divided  into  two  sub-sites  for  analysis.  Sub¬ 
site  1  (see  Figures  10  and  II)  consists  of  an  area  in  the  Thames  River 
Estuary  slightly  northwest  of  Ramsgate,  England.  Sub-site  2,  a  deeper 
water  area,  is  a  portion  of  the  English  Channel  northwest  of  Dunkirk, 
France  (see  Figures  12  and  13)  .  In  sub-site  1,  where  a  93 %  correlation 
was  found,  numerous  sand  and  mud  banks  are  present,  and  the  surface 
features  occurred  over  or  near  these.  In  sub-site  2,  where  a  100% 
correlation  was  found,  the  features  were  present  over  a  series  of 
linear-shaped  shoals  in  the  middle  of  the  English  Channel. 

Figure  14  shows  SAR  imagery  from  over  Study  Site  E,  North  Rona 
Rock.  Here,  the  oval-shaped  surface  feature  corresponds  to  a  shoal 
area  which  surrounds  this  island,  as  can  be  seen  in  Figure  15. 
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gure  11.  bathymetric  Chart  of  Thames  River  Kstuary  (see 
figure  JO  for  corresponding  SAK  imagery). 
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figure  J  2  lor  corresponding  SAH  imagery). 


Figure  14.  SAR  Inage  Surface  Feature  near  North  Rona  Rock, 
Scotland,  from  SEASAT  Rev.  762,  19  August  1978 
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Figure  15.  Combination  of  Figure  14  with  Bathymetric 
Chart  of  North  Rona  Rock.  Area  (depths  in 
fathoms) . 
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PHOTOGRAPHIC  ENHANCEMENT  AND  OPTICAL  PROCESSING  ANALYSIS 

As  mentioned  in  the  previous  section,  some  of  the  bottom-related 
surface  signatures  observed  on  the  SEASAT  imagery  were  quite  subtle  in 
nature,  while  others  were  quite  distinct.  Study  Sites  A  and  B  fell  into 
the  first  category,  and  Sites  C,  D,  and  E,  the  latter.  Because  much  previous 
work  has  been  conducted  at  Study  Site  A  (DMA-Bahamas  Photobathymetric  Calibra¬ 
tion  Area)  in  the  field  of  remote  sensing  bathymetry,  excellent  ground  truth  exists 
for  this  area.  Therefore,  the  decision  was  made  to  both  photographically 
enhance  and  conduct  special  optical  processing  of  the  SEASAT  SAR  data 
from  this  area.  The  purpose  of  these  special  techniques  was  to  see  if 
additional  bottom-related  surface  signatures  could  be  seen  on  the 
imagery,  as  well  as  make  any  present  (but  subtle)  features  more  distinct. 

There  are  several  reasons  for  attempting  special  techniques  in  pro¬ 
cessing  SEASAT  SAR  data.  The  first  concerns  exploiting  the  proper  dynamic 
range  window  for  oceanographic  applications.  Radar  backseat  ter  (a^)  is 
the  wav  in  which  the  radar  return  signals  are  discussed.  These  back- 
scatter  measurements  are  commonly  converted  to  decibel  (dB)  values. 

Radar  backscatter  signals  have  a  dynamic  range  covering  up  to  60  dB  (see 
Figure  IS) .  The  radar  signal  film  is  processed  onto  output  (image) 
film,  which  has  a  dynamic  range  of  only  17-20  dB .  One  can  either  com¬ 
press  the  data  or  only  exploit  a  certain  part  of  the  data,  depending  on 
the  application,  when  recording  Che  data  on  output  film.  Radar  signals 
have  to  be  compressed  before  putting  them  on  film. 

One  can  assign  a  certain  wide  band  (i.e.,  dB  range)  of  data  to  a  narrower 
band,  thus  enabling  the  recording  of  essentially  the  entire  band  (or  dB) 
width  onto  the  recording  medium.  In  either  case,  information  is  lost. 

With  J?L  survey  processed  data,  one  does  not  know  what  part  of  the  entire 
dynamic  range  of  the  SEASAT  data  has  been  displayed  or  what  sort  of  data 
compression  has  been  used.  Custom  optical  processing  of  SAR  signal  film 
allows  exploitation  of  the  proper  dynamic  range  for  oceanographic  studies. 
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Surface  features  which  are  present  on  SAR  imagery  over  ocean  areas 
are  most  likely  due  to  a  physical  oceanic  process.  The  term  "physical 
oceanic  process"  implies  a  motion-related  process.  Previous  SAR  studies 
(see  Shuchman  and  Zelenka,  1978  and  Teleki,  et  al.,  1978)  have  indicated 
that  the  SAR  ocean  surface  imaging  mechanism  is  sensitive  to  motion 
effects  of  the  radar  backscatterers .  Under  certain  circumstances  these 
motion-related  distortions  can  be  corrected  when  processing  the  SAR 
signal  histories  and  additional  surface  features  become  visible.  Mo 
such  corrections  were  made  when  the  JPL  survey  processed  imagery  was 
produced. 

JPL  survey  processed  film  from  SEASAT  Revolution  407  underwent 
special  photographic  processing  in  an  attempt  to  enhance  those  features 
present  and  bring  out  previously  undetected  signatures.  First,  photo¬ 
graphic  prints  were  made  using  low  and  high  contrast  photographic  paper. 
An  example  of  low,  normal  and  high  contrast  images  from  the  area  around 
Bimini  Island  are  presented  in  Figure  17.  Next,  three  to  four  times 
enlargements  were  produced  and  examined  for  possible  bottom-related 
features.  Upon  examination  of  these  sets  of  images,  it  was  concluded 
that  both  these  techniques  provided  little  or  no  more  information  than 
the  survey  processed  imagery.  This  result  was  not  entirely  unexpected. 
The  general  shape  of  the  bottom  in  this  area  resembles  a  flat  or  gently 
rolling  feature.  Therefore,  any  water  flow  over  this  bottom  type 
would  generally  be  less  turbulent  than  say  over  distinct  shoals  (such 
as  around  Nantucket  Island). 

Next,  SAR  signal  film  from  SEASAT  Revolution  651  underwent  special 
optical  processing.  Three  analyses  were  run:  1)  adjustments  were  made 
so  the  optimum  dynamic  range  for  oceanographic  applications  was  dis¬ 
played  on  the  output  film;  2)  different  recording  films  with  different 
Gammas  (contrasts)  were  used;  and  3)  azimuth  focus  corrections  were 
used  to  compensate  for  moving  ocean  surfaces  imaged  by  the  SAR. 

One  test  was  run  where  the  optimum  range  was  exploited  while  test¬ 
ing  various  film  types  which  differed  as  to  their  Gammas.  (Gamma  is  a 
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measure  of  the  contrast  of  the  film,  the  higher  the  Gamma,  the  greater 
the  contrast) .  By  processing  Kodak  EK3414  film  on  an  automatic 
Versamat  Machine,  the  hi^  ast  Gamma  of  1.7  was  obtained.  By  process¬ 
ing  Kodak  panchromatic  on  the  Versamat,  a  Gamma  of  0,8  was  obtained.  3v 
hand  processing  the  EK3414  film,  a  low  Gamma  of  0.6  was  obtained. 

Figure  18  gives  an  example  of  each  of  the  different  Gammas  used. 
Upon  examination  of  these  films,  it  was  again  observed  that  no  addi¬ 
tional  surface  anomalies  could  be  detected. 

A  study  to  see  if  an  azimuth  focus  correction  results  in  enhanced 
imagery  was  conducted  next.  The  Depth  of  Focus  (D.F.)  for  this  imagerv 
(see  Shuchman  and  Zelenka,  1978)  is: 


2  vf 
D.F.  -  f-  (•  £ 


fa.  2 


(1) 


where:  (Illumination  wavelength)  »  0.6328  x  10~6  meters, 


~  =  207,000  ’  where  Vf  is  the  speed  of  the  signal  film, 

SAX 

and  VSAT  is  che  veiocity  of  the  SEASAT  satellite) , 
c  (Azimuth  Resolution)  =  10  meters 

cl 

M  (Processor  Azimuth  Demagnification)  =3.5 

A  depth  of  focus  of  0.6  mm  or  2.4  x  10  d  inches  was  calculated  for 
the  data.  Based  on  this  value,  a  series  of  images  was  made  where  che 
azimuth  focus  shift  was  changed  in  steps  of  approximately  2  x  D.F. 

_3 

(5.0  x  10  inches).  A  series  of  output  films  were  made  using  6  dif¬ 
ferent  focus  settings  +4  D.F.  (+.01  inches),  +2  D.F.  (+.005  inches), 
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c)  Gamma  of  0.6  (low  contrast) 


Figure  18.  Examples  of  SAR  Imagery  Made  Using 
Different  Recording  Film  Types 
During  Optical  Processing  (Rev.  651  Data; 
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0  D.F.  (0.00  inches),  -2  D.F.  (-0.005  inches),  -4  D.F.  (-0.01  inches), 
and  -6  D.F.  (0.015  inches) .  Several  of  these  variable  focus  setting 
films  are  illustrated  in  Figure  19.  Again,  no  additional  surface 
features  were  detected  on  the  specially  focused  data. 
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RELATIVE  RADAR  BACKSCATTER  MEASUREMENTS 

In  order  to  more  rigorously  analyze  the  relationship  between  the 
observed  SAR  surface  features  and  corresponding  bathymetric  information, 
measurements  of  the  radar  backscatter  return  were  compared  to  water 
depth  profiles.  These  measurements  were  obtained  in  two  ways.  Back¬ 
scatter  measurements  were  obtained  by  placing  a  scanning  device  in  the 
output  plane  of  the  optical  processor.  The  output  plane  of  the  optical 
processor  is  the  position  where  the  film  (output)  is  placed  to  record 
the  output  imagery.  Scanning  at  the  output  plane  of  the  processor 
enables  one  to  avoid  the  log  transformation  of  a  film,  and  also  retain 
the  full  dynamic  range  of  an  image.  The  scanning  device  used  in  this 
study  utilized  a  source  aperture  that  corresponded  to  an  ocean  area  of 
350  x  350  M  .  This  technique  requires  the  use  of  SAR  signa.  histories 
of  a  study  area.  The  other  method  for  obtaining  backscatter  measure¬ 
ments  was  through  the  use  of  JPL  digitally  processed  CCT  data.  For 
convenience,  the  backscatter  measurements  were  obtained  in  transect 
form,  which  could  then  be  compared  to  the  bathymetry  of  a  test  site. 
Optical  processor  scans  were  obtained  for  three  of  the  study  sites 
(A,  D,  E)  and  JPL  CCT  data  was  used  on  Study  Site  C. 

It  should  be  noted  that  radar  backscatter  values  presented  in 
this  section  are  relative  radar  backscatter  values.  The  SEASAT  SAR 
was  not  calibrated  in  terms  of  producing  actual  ~ radar  cross  section 
values.  The  relative  radar  backscatter  scans  have  not  been  corrected 
for  affects  of  antenna,  range  power  loss,  sensitivity  time  control 
(SIC),  or  resolution  size.  These  corrections  are  rather  time  con¬ 
suming,  but  should  be  considered  in  future  evaluations. 


43 


^ER! _ 

Figure  20  summarizes  che  relative  radar  backscacter  measurements 
made  in  the  optical  processor  of  Study  Site  A,  located  near  Bimini 
Island.  In  Scan  1,  we  can  see  an  almost  50%  increase  in  relative  radar 
backscatter  when  we  compare  the  deep  water  area  with  the  shallow  water 
area.  In  Scan  2,  we  note  a  gradual  increasing  trend  as  we  pass  from 
deep  to  shallow  water,  but  che  difference  in  backscatter  measurements 
between  deep  and  shallow  water  is  not  as  dramatic  as  with  Scan  1.  The 
ocean  process  responsible  for  the  difference  between  Scan  1  and  2  is 
not  well  understood  at  present.  Speculations  as  to  the  cause  of  the 
differences  between  Scan  1  and  2  will  be  discussed  in  the  next  section. 

Figure  21  summarizes  che  optical  processor  backscatter  scans  made 
for  sub-sice  1  of  Study  Site  D  (the  Thames  River  Estuary  region) ,  and 
Figure  22  summarizes  the  scans  made  for  sub-site  2  (the  English  Channel 
region) .  Both  secs  of  scans  show  a  close  correlation  between  a  surface 
signature  (a  sharp  rise  or  drop  in  backscatter)  and  a  distinct  water 
depth  change.  It  should  be  noted  chat  these  scans  were  made  from  non- 
geomecrically  corrected  imagery,  and  as  such,  only  represent  approximate 
correlations.  In  Figure  21,  a  sharp  increase  in  water  depth  is  associated 
with  an  increase  in  backscacter.  In  Figure  22 »  a  decrease  in  water 
depth  is  associated  with  an  increase  in  backscatter.  In  scan  a  of 
Figure  22  »  a  slight  backscatter  drop  occurs  just  prior  to  the  sharp 
rise  over  the  shoals. 

Figure  23  summarizes  the  optical  processor  backscatter  made  at 
Study  Site  E.  We  have  already  noted  the  close  correlation  between  the 
SAR  observed  surface  features  and  the  bathymetry  of  this  area  (see 
Figure  15').  It  should  be  noted  that  the  very  large  peaks  in  Scans  B 
and  D  of  Figure  23  are  from  North  Rona  Rock,  and  do  not  represent  an 
ocean  surface  variation.  It  is  clear  from  Figure  22,  however,  chat  the 
relative  backscatter  from  the  shoal  area  is  higher  than  from  the  surrounding 
deep  water  areas  (i.e.,  as  the  water  depth  decreases,  the  backscatter 
increases).  The  underlying  oceanographic  processes  responsible  for  the  back¬ 
scatter  variations  (noted  above)  over  the  shoal  regions  are  not  yet  under¬ 
stood  completely. 
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The  JPL  digitally  processed  SAR  .image  from  Nantucket  Island 
(Revolution  880)  is  shown  in  Figure  24.  The  data  has  25  m  resolution 
and  represents  an  area  of  100  x  100  km.  If  Figure  24  (a  positive 
print)  is  compared  to  Figure  8,  a  photographic  negative  of  JPL 
optically  processed  SAR  data,  the  better  quality  of  Figure  24  is 
apparent. 


figure  25  shows  the  location  where  relative  radar  backscatter 
measurements  were  digitally  obtained. 

The  digital  processing  consisted  of  extracting  5  lines  of  data 
at  each  of  the  shown  transects.  These  were  then  averaged  to  form  a 
single  line,  which  in  effect  was  smoothing  the  data  in  the  azimuth 
dimension.  Geometric  corrections  were  then  performed  using  techniques 
outlined  by  Shuchman,  et  al.  (1978).  Appendix  B  presents  a  summary 
of  the  geometric  correction  used  along  with  a  verification  of  the 
algorithm.  Next,  the  data  was  low-pass  filtered  utilizing  ->  variable 
aperture  sliding  window  technique  to  smooth  the  data  in  the  range  direction. 
These  scans  were  then  all  normalized  by  subtracting  from  each  pixel  the  maxi¬ 
mum  value  obtained  and  converting  them  to  dB.  Bathymetry  cross-sections  were 
also  produced  corresponding  geographically  to  the  backscatter  transects.  The 
bathymetry  cross-section  data  were  obtained  from  N0AA  Chart  No.  13237. 

The  digital  backscatter  transects  and  bathymetric  cross-sections 
are  shown  in  Figures  26-23.  Examination  reveals  an  inverse  correla¬ 
tion  of  backscatter  return  and  water  depth,  that  is,  the  larger  the 
depcn,  the  smaller  the  corresponding  backscatter,  and  vice  versa. 

It  can  be  observed  from  Figures  26-23  that  the  correlation  of  radar 
backscatter  appears  to  be  quite  good.  In  general  the  radar  signal 
increases  3  dB  (i.e.,  a  factor  of  two)  over  the  topographic  ridges 
present  on  the  sea  floor  in  the  Nantucket  Shoals  area.  The  profiles 
presented  in  the  figures  indicate  that  the  radar  backscatter 
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FIGURE  -■* .  SAR  Imagerv  maae  from  .'PL  Produced  Digital  Data,  Nantuc.rer 
:iidna,  oEASAt  Pev.  380,  IT  Augusc  1973 
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value  increases  directly  above  the  topographic  ridges,  as  opposed  to 
slightly  preceeding  or  following  the  position  of  the  ridge.  Un¬ 
fortunately,  the  accuracy  limitations  of  the  bathymetric  chart  pre¬ 
vented  the  above  posed  theory  from  being  examined  in  greater  detail. 

Plans  are  currently  being  finalized  to  further  study  the  surface 
features  observed  in  the  Nantucket  Shoals  area.  This  experiment 
(named  SEBEX  for  Surface  Expression  of  Bathymetry  Experiment)  will 
bring  together  oceanographers  as  well  as  SAR  engineers  to  investigate 
both  the  spatial  occurrence  of  these  features  as  well  as  the  SAR/ sea 
surface  interactions  which  cause  them. 
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DISCUSSION 

This  study  has  shown  that,  in  the  limited  number  of  areas  studied, 

SEASAT  synthetic  aperture  radar  (SAR)  image  signatures  can  be  related  to 
bottom  topographic  features.  The  depth  related  surface  signatures  appear  to 
be  present  on  the  SAR  data  only  when  currents  or  gravity  waves  are  present. 
Additionally,  the  orientation  of  the  waves  or  currents  in  respect  to  the 
bottom  feature  appears  to  be  important.  The  effect  of  local  wind  (including 
orientation  and  magnitude)  on  the  visibility  of  bottom  related  signatures 
on  SAR  imagery  is  not  known  at  present  and  needs  to  be  studied  in  detail. 

It  is  known,  however,  that  a  minimum  wind  (on  the  order  of  6  knots)  is 
needed  to  generate  the  capillary  wave  field  responsible  for  radar  back- 
scatter  from  the  sea  surface. 

The  SAR  signatures  are  believed  to  be  an  indirect  result  of  a  hydro- 
dynamic  interaction  between  an  ocean  process  (or  air/ocean  process)  and 
the  bottom  shapes  and  depth  features.  This  concept  can  be  better  under¬ 
stood  if  an  example  is  presented  where  a  distinct  ocean  process  results  in 
an  alteration  of  the  sea  surface,  which  in  turn  can  be  imaged  by  a  satellite 
sensing  system  such  as  SEASAT.  The  illustration  used  will  be  the  internal 
waves  presented  earlier  in  Figure  1.  Figure  29  shows  a  schematic  (after 
LaFond  and  Cox,  1962)  where  part  of  the  energy  field  from  an  internal  wave 
is  being  transmitted  towards  the  surface.  This  energy  in  turn  creates  an 
area  where  the  sea  surface  is  rougher  than  adjacent  areas.  An  area  cf 
increased  roughness  is  caused  by  the  interaction  between  the  internal  waves 
and  the  surface  capillary  waves.  The  internal  wave  energy  at  the  ocean's 
surface  causes  a  "bunching"  of  the  capillary  waves.  This  bunching  of 
the  capillary  waves  creates  a  situation  where  there  are  more  radar  back- 
scatterers  per  unit  area  than  in  adjacent  areas,  thus  the  increase  roughness 
of  this  area.  The  altering  rough-s_.ooch  pattern  on  the  ocean's  surface  (above 
the  internal  wave)  can  then  be  imaged  by  the  SEASAT  SAR,  as  they  were  in 
Figure  1.  This  basic  concept  (that  an  ocean  process  can  result  in  altering 
the  surface  roughness  to  a  degree  which  is  detectable  on  SAR  imagery)  can 
be  expanded  to  include  effects  such  as  currents  and  gravity  wave  motion. 
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FIGURE  29.  Schematic  of  Che  Structure  of  Progressive  Internal  Wave 
Motion  Along  a  Sharp  Thermocline  Illustrating  Alteration 
of  Surface  Roughness  (after  Lafond  and  Cox,  1962) 
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The  principle  in  Imaging  any  ocean  surface  with  a  radar  is  that 
the  backscatter  of  microwave  energy  (echo)  received  by  the  radar  re¬ 
ceiver  contains  information  on  the  roughness  characteristics  (shapes, 
dimensions,  and  orientations)  of  the  reflecting  area.  Parameters  that 
influence  the  radar  return  received  from  ocean  surfaces  include  the  motion  of 
the  scattering  surfaces,  the  so-called  speckle  effect,  system  resolu¬ 
tion  and  non-coherent  integration,  as  well  as  contributions  attributable 
to  wind,  waves,  surface  currents  and  surface  tension. 

The  final  output  of  a  SAR  system  is  usually  in  the  form  of  a  two- 

dimensional  display  of  relative  radar  backscatter  intensity.  The  basic 

backscattering  mechanism  is  believed  to  be  Bragg  scattering  (Wright, 

1966);  that  is,  transmitted  radar  energy  with  wave  number  K  interacts 

in  a  resonant  fashion  with  ocean  surface  waves  with  wave  number  K  such 

w 

that 

K  -  2K  sin  8,  (2) 

w 

where  K  *  2tt/L  and  K  *  2ir/X  are  the  wave  numbers  and  L  and  X  the 
w 

wavelengths,  respectively,  of  the  surface  waves  and  the  radar,  and  8 
is  the  incidence  angle.  For  the  SEASAT  parameters,  the  Bragg  wavelength 
is  between  30  and  40  centimeters.  Shuchman,  et  al.  (1978a)  showed  that  a 
3ragg-Rice  scattering  equacion  satisfactorily  explained  the  radar  backscatter 
return  from  SAR  using  data  collected  during  the  Marineland  experiment  (Shemdin, 
et  al.,  1978).  It  should  be  noted  that  radar  data  from  large  ocean  areas 
(lxl  km)  were  averaged  in  that  analysis.  Thus  (based  on  the  above),  the 
principle  radar  reflectivity  mechanism  of  imaging  ocean  surfaces  is  via 
capillary  and  small  gravity  waves  (Raney  and  Shuchman,  1978).  Any  alteration 
of  the  structure  of  these  waves  will  result  in  a  backscatter  variation 
in  the  recorded  signal  histories. 

In  the  following  discussion,  we  will  present  some  of  the  proposed 
theories  that  attempt  to  explain  the  SAR  signatures  observed  in  each  of 
the  test  sites.  This  section  of  the  report  cannot  be  considered  com¬ 
plete  as  all  proposed  theories  are  not  presented,  and  much  more  work  is 
necessary  in  this  area.  The  SAR  data  cannot  be  effectively  used  as  an 
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indicator  of  bottom  topographic  information  until  the  limitations  of  the 
SAR  data  are  understood.  The  limitations  of  using  SAR  imagery  to  obtain 
topographic  information  by  SAR  will  be  best  understood  by  development 
of  a  theory  of  the  imaging  mechanism. 

The  least  distinctive  (in  terms  of  bottom  bathymetry)  signatures 
occurred  in  the  SAR  imagery  from  Sites  A  and  B,  which  were  located  in 
the  Bahama  Island  group.  The  only  surface  patterns  which  were  believed 
to  be  ba thyme trically  related  are  along  the  edge  of  either  the  Little 
or  Great  Bahama  Banks  where  the  bathymetry  goes  from  shallow  to  very 
deep  water  over  a  short  ground  distance.  In  no  case  was  the  entire 
edge  of  either  of  these  barucs  reflected  on  the  SAR  imagery. 

A  possible  mechanism  for  the  presence  of  SAR  signatures  in  these 
study  sites  lies  in  the  SAR's  ability  to  detect  surface  currents  (see 
Shuchman,  et  al.,  1979a).  An  example  is  the  nearshore  region  near 
Bimini  Island  (see  Figure  8),  where  a  one  to  two  knot  (0.5  to  1.0  m/s) 
northerly  current  exists  as  a  result  of  the  Gulf  Stream.  Because  of  the 
Banks,  this  current  would  not  extend  eastward  due  to  frictional  forces, 
and  this  current  differential  (or  current  shear)  is  then  detected  on  the 
SAR  imagery. 

Only  the  western  portion  of  the  Great  and  Little  Bahama  Banks  have 
been  studied  here.  In  this  area,  very  few  distinct  bottom  features  exist 
(ocher  chan  the  islands  and  shoals  along  the  edges  of  these  bands  and  a 
sharp  drop-off  at  their  edges).  Consider  the  above  observation  with  the 
absence  of  significant  tidal  currents  (probably  less  than  one-half  knot  or 
.25  m/s),  at  che  time  of  che  SEASAT  overpasses,  and  a  situation  exists 
where  there  is  little  possibility  for  an  interaction  between  a  bottom 
feature  and  an  ocean  process  to  occur.  Other  authors  have  reported  the 
presence  of  distincc  bottom  features  on  Landsat  imagery  (Hine,  1977  and 
Defense  Mapping  Agency,  1980),  and  on  SEASAT  imagery  (Sherman,  1980). 

An  example  of  one  such  area  is  presented  in  Figure  30,  which  is  from 
SEASAT  Rev  1411,  collected  on  3  October  1978.  Figure  30  was  collected 
over  the  Tongue  of  the  Ocean,  a  prominent  deep  water  trench  which  extends 
into  the  Little  Bahama  Bank.  As  can  be  seen,  prominent  surface  features 
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exist,  which  are  believed  to  be  bottom  related.  SEASAT  made  at  least  five 
passes  over  this  area  and  this  represents  a  study  site  for  future  studies. 

As  noted  earlier  from  Figures  21  and  22,  we  have  two  opposite 
situations  occuring  at  Site  D.  In  the  Thames  Estuary  Region  (Figure  21) 
an  increase  in  water  depth  is  associated  with  an  increase  in  backscatter. 

In  the  English  Channel  Region,  a  decrease  in  water  depth  is  associated 
with  an  increase  in  backscatter.  One  would  expect  that  two  separate 
oceanographic/bottom  interactions  are  responsible. 

When  interpreting  SAR  ocean  imagery,  one  associates  a  decrease 
in  radar  backscatter  with  decreased  water-surface  roughness  and  vice- 
versa.  The  Thames  River  Estuary  is  an  area  of  very  shallow  (<lm) 
tidal  flats  and  it  is  quite  possible  that  the  capillary  waves  over 
these  flats  are  sufficiently  suppressed  to  result  in  a  significant 
decrease  in  surface  roughness. 

The  backscatter  variations  observed  in  the  English  Channel  sub-site  are 
oeiieved  to  be  tidal  current  induced.  The  exact  mechanics  of  these  tidal 

overfall  areas  will  be  discussed  in  greater  detail  below.  Figure  31 
shows  the  tidal  currents  which  were  present  during  the  SEASAT  Rev.  762. 

Note  that  1  to  2.3  knot  (0.5  to  1.25  m/s)  currents  existed.  Unfortunately, 
the  wind  conditions  on  the  day  of  the  overpass  are  not  known  at  the 
time  of  this  writing.  Figure  32  is  from  a  Landsat  study  of  the  Test 
Site  D  performed  by  Viollier  and  Baussart  (1979) .  Of  particular  interest 
to  this  study  is  the  MSS  Band  7  (800  to  llOOnm)  imagery  in  which 
there  is  essentially  no  penetration  into  the  water  column  and  hence 
any  patterns  over  water  are  the  result  of  the  sun's  reflection  off  the 
capillary  waves  on  the  water  surface.  In  the  12  June  1975  Band  7 
imagery  (on  Figure  32),  the  same  linear  patterns  are  present  in  the 
center  of  the  English  Channel  as  were  present  in  the  SEASAT  imagery 
(Figure  12) of  the  same  area.  Tidal  current  conditions  were  essentially 
the  same  during  the  collection  of  the  12  June  1975  Landsat  and  the  Rev. 

762  SEASAT  data,  with  a  strong  tidal  flow  to  the  southwest.  The  30 
March  1977  Landsat  imagery  in  Figure  32  was  collected  at  a  time  when 
the  tidal  flow  was  opposite  in  direction  (towards  the  northeast)  to 
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lanasac  Imagery  or  Study  Sica  0,  Suo-site  ’  (English  Channel) 
from  Vioiiier  ana  3aussart  [1979!  vach  oermission  of  the  authors) 


Che  Cides  present  on  12  June  1975.  No  depth  related  signatures  can 
be  observed  on  the  30  March  1977  Landsat  data.  It  should  be  noted 
Chat  very  little  wind  was  present  during  the  March  imagery  and  the  sun 
illumination  was  low,  a  poor  set  of  conditions  for  sun  glitter  off  the 
water. 

The  Landsat  imagery  from  this  study  is  important  from  several 
standpoints.  It  indicates  that  the  depth  related  features  observed 
on  the  SEASAT  imagery  were  not  just  coincidental  to  the  time  of  the 
SEASAT  overflight.  Also,  the  comparison  of  Landsat  imagery  from  two 
dates  when  opposite  tidal  flows  and  different  wind  conditions  were 
present  shows  the  necessity  of  a  distinct  oceanic  physical  process  in 
order  that  the  bottom  related  surface  feature  be  present  on  the  imagery. 

The  possible  mechanism  for  the  surface  features  in  the  shoal  areas 
around  North  Rona  Rock  (Study  Site  E,  see  Figures  14  and  15)  is  not  yet 
clear,  but  their  occurrence  directly  over  a  shoal  area  which  is  surrounded 
by  deeper  water  leaves  very  little  doubt  that  they  are  related  to  the 
bottom  bathymetry  of  the  area. 

One  possible  mechanism  for  the  presence  of  surface  anomalies  on 
SAR  imagery  over  North  Rona  Rock  is  a  non-linear  interaction  between 
gravity  waves  and  surface  capillary  waves  due  to  the  sharp  rise  (^40 
meters)  in  the  bottom  in  the  shoal  region.  This  sharp  rise  causes  a 
drastic  alteration  in  the  gravity  wave  structure  which  results  in  both 
wave  refraction  and  a  gravity-capillary  wave  interaction.  The  surface 
signatures  might  possibly  be  tidal-current  induced,  although  the  tides 
in  this  area  at  the  time  of  the  SEASAT  overpass  were  probably  less  than 
one-half  knot  (.25  m/s)  in  an  easterly  direction. 

Figure  33  summarizes  the  tidal  and  wide  conditions  present  at  the 
time  of  the  SEASAT  Rev.  880  over  Nantucket  Island  (Study  Site  C) . 

Weather  records  were  consulted  to  see  if  any  major  climatic  events 
occurred  on  or  before  the  SEASAT  overpass.  The  region  was  dominated 


by  a  high  pressure  center  during  and  several  days  prior  to  the  SEASAT 
overpass.  Wind  speeds  were  generally  light,  in  the  range  of  5-10  knots 
(2.5  to  5.0  m/s)  from  the  north. 

It  is  believed  that  the  surface  patterns  in  the  inshore  regions 
of  Rev.  880  are  the  result  of  an  interaction  between  the  tidal  currents 
and  the  shoals  in  the  Nantucket  Island/Cape  Cod  region.  This  notion 
is  again  supported  by  observations  made  from  Landsat  imagery.  Figure 
34  is  a  Band  6  (.7-. 8  um)  image  collected  on  17  July  1974.  Like  Band 
7,  Band  6  has  very  little  water  penetration,  and  therefore,  water 
features  present  on  the  imagery  are  a  result  of  reflected  sunlight  from 
the  capillary  waves  on  the  ocean  surface. 

The  relative  radar  backscatter  scans  (Figures  26-28)  made  for  this 
area  showed  the  close  correlation  between  backscatter  and  water  depth 
in  this  area.  The  scans  could  not  positively  answer  the  question  whether 
the  increase  in  radar  backscatter  (associated  with  a  decrease  in  water 
depth)  occurred  over  the  topographic  "bump"  in  bottom  or  slightly  pre- 
ceeding  or  after  the  features.  It  is  assumed  that  the  large  rise  in  radar 
backscatter  observed  on  the  Nantucket  Shoals  data  occurs  slightly  after  the 
ridge  in  bottom.  The  oceanographic  cause  of  the  increase  in  backscatter  is 
commonly  referred  to  as  a  tidal  overfall  and  has  been  well  documented  by 
Stewart  and  Jordan  (1964).  Stewart  and  Jordan  studied  the  underwater  and 
ridges  of  Georges  Shoal  which  is  located  approximately  120  km  east  of  the  Nan¬ 
tucket  test  area.  The  Georges  Shoal  is  very  similar  to  the  Nantucket  Shoals 
in  that  an  elongated  bank  of  shallow  water  (over  the  shoals)  lies  next  to  an 
area  of  deeper  water  (usually  10  to  15  meters  deeper).  Stewart  and  Jordan 
reported  that  when  a  strong  (in  excess  of  1-1/2  kts  or  Im/s)  tidal  current 
flowed  normal  (perpendicular)  to  the  ridges  on  Georges  Shoal,  an  overfall  is 
created  as  the  large  volume  of  fast  moving  water  is  forced  upward  and  over 
the  ridge.  An  actual  rise  of  the  sea  surface  was  observed  by  Stewart 
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and  Jordan  just  down  current  from  the  crest  of  the  underwater  ridge. 
The  overfall  "wave"  was  stationary  relative  to  the  bottom  feature 
which  was  creating  it  (along  with  the  tidal  current) .  These  overfall 
waves  were  only  present  during  times  when  the  tidal  currents  were 
strong  and  normal  to  the  ridges .  The  overall  waves  are  short-period 
waves  and,  typically,  a  short  period  breaker  line  can  develop  (Stewart 
and  Jordan,  1964) . 

From  Figure  33 ,  it  is  clear  that  the  tidal  currents  are  generally 
perpendicular  to  the  underwater  topographic  features.  Therefore  a 
similar  situation  to  that  reported  by  Stewart  and  Jordan  (1964)  could 
have  occurred  at  this  study  site,  which  would  account  for  the  large 
scale  surface  features  on  the  SAR  imagery. 

The  observations  made  from  the  Landsat  imagery  in  Figure  34  were 
also  consistent  with  Stewart  and  Jordan's  work.  The  tides  on  the  date 
the  Landsat  data  was  collected  were  similar  in  magnitude  and  direction 
as  at  the  time  of  the  SEASAT  overpass  (Reed,  1980) . 
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CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  helped  to  demonstrate  that  synthetic  aperture 
radar  (SAR)  data  has  the  potential  to  detect  bathymetric  features  in 
shallow  water  regions.  The  detection  process  is  via  a  hydrodynamic 
interaction  between  currents  and/or  waves  with  the  bottom  topographic 
feature  which  alters  the  surface  roughness  of  the  ocean.  The  SAR  is 
very  sensitive  to  surface  roughness  changes  of  the  ocean  surface,  and 
because  microwave  sensors ,  like  SARs ,  do  not  penetrate  more  than  a 
few  centimeters  into  the  ocean,  the  depth-related  signatures  are  ocean 
surface  effects. 

The  greatest  potential  exists  for  SAR  detection  of  bottom  features 
in  areas  where  there  are  shoals  with  currents  flowing  over  them.  Pre¬ 
sently,  the  SAR  data  can  only  be  used  as  an  indicator  of  shoal  areas, 
but  perhaps  this  information  can  be  used  to  ascertain  whether  a 
bathymetric  chart  needs  updating  or  revision.  This  could  be  done  by 
correlating  the  SAR  signatures  to  a  digitized  bathymetric  chart. 

Whether  a  theory  and  a  resulting  processing  algorithm  can  be  developed 
that  relates  the  radar  backscatter  signatures  observed  in  this  report 
directly  to  actual  water  depth  is  not  known,  and  should  be  an  area  of 
future  research. 

Areas  where  distinct  bathymetric  features  occur  in  coastal  areas 
and  in  which  physical  oceanic  processes  occur  (strong  current  or  wave 
action)  present  the  best  potential  for  bathymetry/ SAR  correlations. 

In  the  present  study,  the  Nantucket  Shoals,  English  Channel,  Thames 
River  Estuary  and  North  Rona  Rock  shoals  represented  such  areas. 
Although  portions  of  the  edge  of  the  3ahama  Banks  were  discernible 
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on  che  SAR  imagery,  bathymetric  information  on  the  Banks  themselves 
was  not;  this  was  most  likely  due  to  the  lack  of  a  physical  oceanic 
process  or  a  distinct  bathymetric  feature.  Recalling  that  the  banks 
are  relatively  flat,  even  if  currents  or  waves  were  present,  they 
need  to  interact  with  a  batVmetric  feature,  which  is  not  present. 
Special  photographic  enhance  nt  and  optical  processing  techniques 
did  not  seem  to  add  any  information  in  the  Bahama  Banks  area. 

Digitally  processed  SAR  data  offers  the  greatest  potential  in 
further  research  because  this  form  of  data  lends  itself  to  computer 
processing  techniques.  Computer  processing  is  necessary  to  geo¬ 
metrically  and  radiometrically  correct  the  large  volume  of  information 
present  in  SAR  data  sets.  Geometric  correction  of  the  data  is  necessary 
for  accurate  comparison  with  localized  bathymetry  data. 

Many  anomalous  surface  features  occur  on  SEASAT  imagery  collected 
over  deep  ocean  areas.  Care  -ust  be  taken  when  interpreting  the 
significance  of  these  features  in  relation  to  bathymetric  information. 
Extensive  study  of  physical  oceanic  processes  and  their  SAR  signals 
have  yet  to  be  conducted.  Investigators  of  coastal  area  remote  sensing 
imagery  must  be  cautious  in  interpreting  SAR  imagery  in  order  that  the 
data  are  correctly  analyzed. 

This  study  consisted  of  a  limited  data  set  in  that  only  a  few 
oceanographic  areas  were  analyzed.  Additionally,  all  the  various 
techniques  available  to  study  the  data  were  not  used  at  each  test 
site.  Further  research  needs  to  be  explored  to  more  fully  exploit 
the  use  of  SAR  data  for  the  detection  of  shallow  water  bathymetric 
features.  They  include: 

a)  examination  of  other  SEASAT  SAR  passes  of  coastal  regions  to 
observe  topographic  bottom  features; 
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b)  examination  of  other  SEASAT  SAR  passes  of  the  five  test 
sites  used  in  this  study.  Hopefully  SAR  data  was  collected 
over  the  five  sites  during  different  environmental  conditions 
(i.e.,  different  wind,  wave  and  current  conditions).  These 
additional  data  sets  could  be  used  to  better  define  the 
limitations  of  SAR  sensing  of  topographic  features. 

c)  utilization  of  more  sophisticated  analysis  techniques  on  test 
areas  that  show  particular  promise. 

Specifically  further  research  in  the  use  of  SEASAT  SAR  data  to 
detect  shallow-water  features  should  include: 

a)  the  quantification  of  environmental  parameters  (such  as  wind, 
•  wave  and  tide  information) ; 

b)  geometric  and  radiometric  (power  loss,  antenna  patterns, 
sensitivity  time  control)  corrections  of  the  SAR  data; 

c)  computer  merging  of  the  SAR  data  with  bathymetric  data  to 
statistically  correlate  the  relationship  between  the  two; 

d)  merging  of  SEASAT  data  with  Landsat  data; 

e)  merging  and/or  comparing  SEASAT  data  from  the  same  area  from 
different  dates; 

f)  based  on  the  above,  explore  the  development  of  theories  and 
algorithms  to  define  the  relationship  between  SAR  detected 
surface  features  and  bottom  topographic  information. 
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BATHYMETRIC  CORRELATIONS 


APPENDIX  A 


BATHYMETRIC  CORRELATIONS 

The  first  step  in  determining  if  SAR  data  could  be  used  to  pro¬ 
vide  bathymetric  information  was  to  determine  whether  or  not  the 
surface  features  observed  on  the  SAR  imagery  occurred  over  or  near 
distinct  bathymetric  features.  This  analysis  was  done  through  a  pro¬ 
cedure  (termed  bathymetric  correlation)  where  a  clear  cell  reproduction 
of  a  hydrographic  chart  (enlarged  or  reduced  to  the  same  map  scale)  was 
overlayed  onto  the  SAR  imagery.  The  water  depth  under  each  feature  was 
then  obtained  from  the  chart,  and  this  sub-feature  depth  was  compared 
to  the  adjacent  area  (within  a  few  hundred  meters)  water  depth.  It 
could  then  be  determined  if  the  feature  on  the  SAR  imagery  occurred 
over  or  near  an  area  where  the  water  was  shoaling  or  getting  deeper. 

For  convenience,  each  Study  Site  was  divided  into  a  group  of 
study  areas.  A  study  area  is  defined  as  the  position  on  the  SAR 
imagery  where  a  feature  believed  to  be  bathvmetrically  related  existed. 
After  all  these  areas  at  a  given  Study  Site  had  been  identified,  the 
clear  cell  reproduction  of  the  chart  of  the  Study  Site  was  then  over¬ 
laid  onto  the  SAR  imagery.  The  next  step  was  to  calculate  the  average 
water  depth  below  each  area  marked  on  the  SAR  imagery.  Then,  the 
average  depth  of  the  water  immediately  adjacent  to  each  study  area  was 
calculated.  If  the  difference  between  these  two  averages  were  greater 
than  50  percent  of  the  shallower  depth,  then  a  distinct  bottom  topo¬ 
graphic  feature  was  said  to  exist. 

The  bathymetric  correlations  for  each  Study  Site  (except  Sices 
3  [Southern  3ahamas ]  and  E  [North  Rona  Rock])  are  presented  in  Tables 


A-i 


A1  to  A5 .  Because  the  number  of  study  areas  at  Sites  B  and  E  were 
limited,  no  tables  were  constructed.  Given  in  each  table  in  this 
appendix  are  the  Study  Area,  the  Average  Water  Depth  for  both  the 
Study  Area  and  the  Adjacent  Area,  and  finally,  a  brief  description 
from  the  SAR  imagery  of  the  characteristics  of  the  study  area  versus 
the  adjacent  area. 


TABLE  A 2  .  BATHYMETRIC  CORRELATION  FOR  SITE  A,  DMA- BAHAMAS  PHOTOBATHYMETR 1 C 
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taule  ac>.  bathymetric  correlation  cor  site  n,  sub-site  2,  English  channel,  seasat  rev.  762 


C 

"3 

>> 

>s 

>i 

w 

2 

-3 

-3 

*3 

*3 

3 

C 

3 

3 

3 

3 

3 

AJ 

aj 

aj 

3 

aj  0) 

4-1 

a; 

AJ  0) 

AJ 

3 

AJ 

aj 

-3 

50 

u 

50  U 

50 

u 

3  U 

X 

u 

2 

"3 

3 

3 

3 

3 

3 

3 

La 

C 

c 

JC 

aj 

^3  aj 

aj 

_3  aj 

JC 

AJ 

3 

CO 

0 

AJ 

X 

•u  X 

AJ 

X 

AJ  X 

*J 

X 

-D 

-a: 

aj 

•pA 

3 

—  cu 

•pJ 

3 

■pj  3 

*— « 

3 

3 

u 

3 

AJ 

3  ^ 

3 

AJ 

3  AJ 

3 

AJ 

C. 

CO 

CO 

u 

z 

MH 

-3 

3 

3 

a 

"3 

•  ro 

** 

*3 

« -0 

*» 

•3 

c 

CO 

u 

L- 

aj 

3 

c  3 

C 

3 

C  3 

c 

3 

50 

< 

AJ 

CO 

CO 

p* 

u 

p«4 

v-  ^ 

u 

fp 

u  — * 

u 

—A 

r-A 

3a 

CO 

•n 

30 

0) 

3  c. 

0) 

3a 

3  a 

3 

o. 

v 

aj 

1A 

AJ 

cu 

aj  C— 

aj 

Q. 

AJ  Cm 

AJ 

50a 

< 

< 

S3 

3 

c 

P-A 

aj 

AJ  iH 

aj 

•pi 

mJ  -pA 

AJ 

•fj 

Q 

3 

3 

u 

3  U 

3 

u 

3  U 

3 

u 

r 

e- 

CO 

1-» 

u 

50 

C* 

a 

Cm 

Cm 

a. 

0 

Z 

2 

u 

CO 

T* 

3 

3 

3 

3 

3 

Aj 

“3 

TO 

*3 

•3 

-3 

CJ 

X 

-3 

"3 

3 

00 

0)  00 

a; 

00 

3  O0 

3 

00 

< 

< 

CO 

CO 

2 

<■■> 

3.  e 

Cm 

c 

CL.  S 

a. 

s 

aa 

—i 

3 

3 

3 

1-4 

— H 

1-4  *w 

*-4 

■W  mJ 

»p^ 

N-» 

as 

u 

3 

jC 

u 

> 

w  > 

> 

u  > 

u 

> 

< 

< 

CO 

3 

_  ~ 

«J 

3 

AJ  3 

AJ 

3 

AJ  CO 

AJ 

3 

*j 

“3 

50 

50  -3 

50 

JZ 

X  -3 

X 

JC 

> 

>N 

3 

C 

a; 

3 

S-' 

u 

C 

c 

3 

3 

a;  3 

3 

3 

3  3 

3 

! 

3 

2 

2 

u 

g 

3 

s  3 

£ 

3 

£  3 

£ 

3 

3. 

AJ 

3 

3 

u 

3  u 

3 

U 

3  U 

3 

u 

cn 

CO 

3 

;o 

CO 

3 

CO  3 

CO 

3 

CO  3 

CO 

3 

S  dJ 

X 

X 

X 

X 

X 

X 

H  CJ 

u 

u 

u 

u 

u 

^a 

a-  z 

3 

3 

3 

3 

3 

3 

aa  a: 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

/— >  myt 

3 

3 

3 

3 

3 

3 

“  aa 

C&a 

£ 

e 

£ 

£ 

£ 

£ 

aa  a- 

ro 

—M 

<r 

r-- 

<  E 

3 

!  ro 

r*' 

CNi 

CN 

CO 

-a* 

<C 

aa 

X 

X 

X 

X 

X 

a: 

X 

u 

U 

u 

u 

U 

< 

u 

3 

3 

V 

3 

3 

3 

AJ 

AJ 

AJ 

AJ 

AJ 

c- 

z 

AJ 

3 

3 

s 

3 

£ 

2 

3 

H 

3 

£ 

CJ 

00 

O 

O 

in 

< 

O 

m 

■* ; 

< 

X 

X 

X 

X 

csa 

X 

X 

La 

La 

La 

u 

u 

u 

3 

3 

3 

3 

< 

3 

3 

w 

j 

«j 

aJ 

«J 

•J 

3 

3 

3 

3 

V 

2 

£ 

s 

5 

rs. 

vO 

00 

A_ 

rvi 

'lo 

<*  t 

O  I 


A- 12 


Terim 


APPENDIX  B 


GEOMETRIC  CORRECTION  OF  SEASAT  SAR  DATA 


APPENDIX  B 


ft 


GEOMETRIC  CORRECTION  OF  DIGITAL  SEASAT  SAR  DATA 


In  order  to  assess  how  closely  the  relative  radar  backscatter 
measurements  present  on  the  JPL  digital  tapes  correlate  to  bathymetric 
information,  it  was  first  necessary  to  geometrically  correct  the  JPL 
data  in  the  range  direction.  The  backscatter  scans  were  made  by 
averaging  five  lines  of  pixels  across  the  range  direction.  The  back¬ 
scatter  scans  were  then  geometrically  corrected  using  the  algorithm 
of  Shuchman,  et  a.,  1978.  In  this  algorithm,  GR  (ground  range  from 
nadir)  is  expressed  as: 


G 


R 


a  cos 


<Pe  +  a)2  -  SR2 

2  O  (pa  +  a) 
e  e 


where : 


(1) 


o  »  radius  curvature  of  the  earth 
e 

a  *  altitude  of  spacecraft 
S  *  slant  range 

This  algorithm  was  tested  using  a  portion  of  SEASAT  Revolution  880 
(see  Figure  1) .  A  transect  was  taken  through  Nantucket  Island  through 
the  city  of  Nantucket.  By  making  a  grey-raap  of  Nantucket  Island,  it  was 
possible  to  locate  where  the  transect  crossed  the  two  shorelines,  and 
the  pixel  numbers  of  these  points  were  recorded.  From  these  pixel 
numbers,  it  was  then  possible  to  locate  their  respective  ground  range 
distances  on  a  computer  printout  from  a  program  based  on  Equation  (1). 
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By  subtracting  the  two  ground  range  distances,  an  estimate  of  the 
transect  distance  over  land  was  produced.  This  estimate  was  17,005 
meters.  The  transect  was  next  located  on  NOAA  Chart  <>13237  (Nantucket 
Sound  and  Approaches) .  This  estimated  distance  for  the  transect  on 
the  chart  was  17,040  meters.  The  difference  between  the  two  estimates 
was  .2%.  Therefore,  the  algorithm  was  judged  to  be  accurate  enough 
for  the  purposes  of  backscatter-bathymetry  comparisons. 
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